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1. Large Transverse Momentum 7t production — from ISR to RHIC 

PHENIX has presented measurements of 71*' production at mid-rapidity in p-p collisions at 
two values of cm. energy ^=200 GeV and 62.4 GeV (Fig. [T]). Some of my younger colleagues 
are amazed at the excellent agreement of Next to Leading Order (NLO) and Next to Leading Log 
(NLL) perturbative Quantum Chromodynamics (pQCD) calculations ^] with the measurements. 
However, this comes as no surprise to me because hard scattering in p-p collisions was discovered 
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Figure 1: (left) PHENIX measurement of invariant cross section, Ed^o/d^p, as a function of transverse 
momentum pr for tt" production at mid-rapidity in p-p collisions at cm. energy = 200 GeV (right) 
PHENIX measurement of 7t^ in p-p collisions at y/s = 62.4 GeV 

at the CERN ISR by the observation of a very large flux of high transverse momentum tt" with a 
power-law tail which varied systematically with the cm. energy of the collision. This observation 
in 1972 proved that the partons of deeply inelastic scattering were strongly interacting. Further ISR 
measurements utilizing inclusive single or pairs of hadrons established that high transverse momen- 
tum particles are produced from states with two roughly back-to-back jets which are the result of 
scattering of point-like constituents of the protons as described by QCD, which was developed 
during the course of these measurements. 



1.1 ISR Data, Notably CCR 1972-73 

The Cern Columbia Rockefeller (CCR) Collaboration (and also the Saclay Strasbourg [Q] 
and British Scandinavian [§] collaborations) measured high pr pion production at the CERN- 
ISR (Fig. HI). The g^^'''^ breaks to a power law at high pr with characteristic dependence. ^ 



' The clear break of the exponential to a power law at ^/s = 200 GeV is shown in the inset of Fig. [l]-(left). 
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The large rate indicates that partons interact strongly EM) with each other, but, to quote 
the authors [^: "Indeed, the possibility of a break in the steep exponential slope observed at 
low pt was anticipated by Berman, Bjorken and Kogut However, the electromagnetic form 
they predict, pj'^F {pj / \/s), is not observed in our experiment. On the other hand, a constituent 
exchange model proposed by Blankenbacler, Brodsky and Gunion [^, and extended by others, 
does give an excellent account of the data." The data fit pj"F{pT /\'^), with « ~ 8. 




Figure 2: (left) CCR ^ transverse momentum dependence of the invariant cross section at five center of 
mass energies, (right) The above data multiplied by p'j, using the best fit value of n = 8.24 ±0.05, with 
F = Ae-^'T, plotted vs pj (= xj /2). 



1.2 Constituent Interchange Model (CIM) 1972 

Inspired by the dramatic features of pion inclusive reactions revealed by "the recent measure- 
ments at CERN ISR of single-particle inclusive scattering at 90° and large transverse momentum", 
Blankenbecler, Brodsky and Gunion [0] proposed a new general scaling form: 

Ep^ = l,F(£L) (11) 

dp^ p} ^/s 

where n gives the form of the force-law between constituents. For QED or Vector Gluon exchange, 
n = 4, but perhaps more importantly, BBG predict n=8 for the case of quark-meson scattering by 
the exchange of a quark (CIM) as apparently observed. 

1.3 First prediction using 'QCD' 1975— WRONG! 

R. F. Cahalan, K. A. Geer, J. Kogut and Leonard Susskind [^] generalized, in their own words: 
"The naive, pointlike parton model of Berman, Bjorken and Kogut to scale-invariant and asymptot- 
ically free field theories. The asymptotically free field generalization is studied in detail. Although 
such theories contain vector fields, single vector-gluon exchange contributes insignificantly to 
wide-angle hadronic collisions. This follows from (1) the smallness of the invariant charge at 
small distances and (2) the breakdown of naive scaling in these theories. These effects should 
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explain the apparent absence of vector exchange in inclusive and exclusive hadronic collisions at 
large momentum transfers observed at Fermilab and at the CERN ISR."^ 

Nobody's perfect, they get one thing right! They introduce the "effective index" ngs{xT,\^) 
to account for 'scale breaking': 

d^a 1 1 
^TTT = 1 rT^(^r) = ^ Gixr) , (1.2) 

where xj = 2pt/\/s. 

1.4 CCOR 1978— Higher pt>7 GeV/c— «eff(^r, V^) ^ 5 = 4++. QCD works! 

The CCOR measurement |^ (Fig. ^ with a larger apparatus and much increased integrated 




Figure 3: a) (left) CCOR transverse momentum dependence of the invariant cross section for 
p + p n'^ +X at three center of mass energies. Cross sections are offset by the factors noted. Open 
points and dashed fit are from a previous experiment, CCRS |]lO[|. b) (right)-(top) Same CCOR invariant 
cross sections plotted vs xj — 2pj j ^fs on a log-log scale, c) (right)-(bottom) neff(x7', y^) derived from the 
combinations indicated. The systematic normalization error at %/s — 30.6 GeV has been added in quadrature. 
There is an additional common systematic error of ±0.33 in n. 

luminosity extended their previous tt" measurement [^, to much higher pj. The pj^ scaling-fit 
which worked at lower pj extrapolated below the higher pj measurements for > 30.7 GeV and 

^There is an acknowledgement in this paper which is worthy of note:"Two of us (J. K. and L. S.) also thank 
S. Brodsky for emphasizing to us repeatedly that the present data on wide-angle hadron scattering show no evidence 
for vector exchange." 
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PT>1 GeV/c (Fig. gi). A fit to tlie new data [|] for 7.5 < /Jr < 14.0 GeV/c, 53. 1 < < 62.4 GeV 
gave Ed^a/dp^ ~ pj^'^'^^'^{l — X7')^^-^^°'^, (including all systematic errors). 

The effective index nef^{xT,^/s) was also extracted point-by-point from the data as shown in 
Fig. |3|b where the CCOR data of Fig. ^a for the 3 values of ^/s are plotted vs xj- on a log-log scale. 
'^eff(^r7\/^) is determined for any 2 values of \/shy taking the ratio as a function of xt as shown 
in Fig. ||c. «eff(-^r, \/^) clearly varies with both -y/i and xt, it is not a constant. For ^/s = 53.1 
and 62.4 GeV, ?ieff(-Xr, \/^) varies from ~ 8 at low xt to ~ 5 at high xj. An important feature 
of the scaling analysis (Eq. pr2[ ) relevant to determining «eff(-^r, \/^) is that the absolute pr scale 
uncertainty cancels! 




Figure 4: (left)-(top) Invariant cross section for inclusive K for several ISR experiments, compiled by 
ABCS Collaboration [[u]]; (left)-(bottom) neffC-^r, V^) from ABCS 52.7, 62.4 GeV data only. There is an 
additional common systematic error of ±0.7 in n. (right)-a) V^(GeV)''-^** x Ed^a/dp^ as a function of 
XT =2/7r/^forthe PHENIX 62.4 and 200 GeV tt" data from Fig. |i]; (right)-b) point-by-point «eff (xr, V^)- 



The effect of the absoulte scale uncertainty, which is the main systematic error in these exper- 
iments, can be gauged from Fig. ^(left)-(top) | [TT| ] which shows the tt" cross sections from several 
experiments. The absolute cross sections disagree by factors of ^ 3 for different experiments but 
the values of ?ieff(-^7-,\/^) for the CCOR [|] (Fig. ||-(right)-(bottom)) and ABCS [jlT]] experiment 
(Fig. ^-(left)-(bottom)) are in excellent agreement due to the cancellation of the error in the ab- 
solute pt scale. The xt scaling of the PHENIX p-p 71° data at = 200 and 62.4 GeV from 
Fig. [I] with nf^fi{xT ,y/s) 6.38 is shown in Fig. ^(right). The log-log plot emphasizes the pure 
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power-law pj dependence of the invariant cross section, Ed^a/dp^ ~ pj" for pj >4 GeV/c, with 
n = 8 . 11 ± 0.05 at = 200 GeV. [g 

The first modern QCD calculations and predictions for high pr single particle inclusive cross 
sections, including non-scaling and initial state radiation were done in 1978, by Jeff Owens and 
collaborators. [12] Jets in 471 Calorimeters at ISR energies or lower are invisible below \/f ~ 



Et < 25 GeV [13]; but there were many false claims which led to skepticism about jets in hadron 



collisions, particularly in the USA. [[T^] A 'phase change' in belief-in-Jets was produced by one 



UA2 event at the 1982 ICHEP in Paris [15], but that's another story. [16] 



2. The major discovery at RHIC-:;r*^ suppression in A+A collisions. 

The discovery, at RHIC, that are suppressed by roughly a factor of 5 compared to point-like 
scaling of hard-scattering in central Au-i-Au collisions is arguably the major discovery in Relativis- 
tic Heavy Ion Physics. In Fig. ^-(left), the data for tt" and non-identified charged particles (/z^) 
are presented as the Nuclear Modification Factor, Raa{pt), the ratio of the yield of tt" (or /j^) per 
central Au-i-Au collision (upper 10%-ile of observed multiplicity) to the point-like-scaled p-p cross 
section: 

fNlJdpTdyNA^ 
""^(^^^ - {T^)d^a;jdprdy ' ^'"^^ 

where (7aa) is the overlap integral of the nuclear thickness functions. The tt" data at nucleon- 
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Figure 5: (left) Nuclear modification factor ^^^(pj) for tt" and in central Au+Au collisions at ^/snn = 
200 GeV [|l7|; (right) Raa{pt) for in Cu+Cu central collisions at ^/s^ = 200, 62.4 and 22.4 GeV [p 
together with theory curves ^ 



nucleon cm. energy ^/smi = 200 GeV are consistent with a constant Raa ~ 0.2 over the range 
4 < < 20 GeV/c, while the suppression of non-identified charged hadrons and tt" are different 
for 2 < < 6 GeV/c and come together for pj > 6 GeV/c. 

A new PHENIX result [18] (Fig. ^(right)) nicely illustrates that parton suppression begins 
somewhere between ^snn=22A and 62.4 GeV for Cu-i-Cu central collisions. This confirms that 
(jet) suppression is unique at RHIC energies and occurs at both ^/s^n = 200 and 62.4 GeV. The 
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suppression is attributed to energy-loss of the outgoing hard-scattered color-charged partons due 
to interactions in the presumably deconfined and thus color-charged medium produced in Au-i-Au 
(and Cu-i-Cu) collisions at RHIC ||^. 

2.0.1 Precise and accurate reference spectra are crucial 

It is important to note that PHENIX did not measure the reference p-p spectrum at ^/s = 
22.4 GeV but used a QCD-based fit [^] to the world's data on charged and neutral pions which was 
checked against PHENIX p-p measurements at 62.4 and 200 GeV using xt scaling (Fig. ^ [jl^]. 
A key issue in this fit is that the data at ^/s = 22.4 GeV were consistent with each other and with 
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Figure 6: Plots for ^^5X^=22. 4, 62.4 and 200 GeV of: a) measured invariant tt" yields in central Cu+Cu 
collisions; b) measured invariant cross sections in p-p collisions at 62.4 and 200 GeV and fit at 
22.4 GeV [|l|]; c) the p-p data and fit from (b) plotted in the form y^iGeV)""" x Ed^a/dp^ to exhibit 
XT scaling with n^n 6.1-6.4, consistent with Fig. 0-(right)-b). 



pQCD [ ]21[ ] except for one outlier which was excluded based on the experience from a previous 
global fit p2| ] to the world's data at = 62.4 GeV where there are large disagreements (recall 
Fig. ^(left)-(top)). The PHENIX measurement of the p-p reference spectrum at 62.4 GeV [Q] 
agreed with the measurements shown in Figs. ^ and ^(left) to within the systematic error of the 



absolute pr scales, but disagreed significantly with the global fit at 62.4 GeV [22] which did not 
attempt to eliminate outliers and which had no basis for adjusting the absolute pr scales of the var- 
ious measurements. In Fig. ^(left), the Raa{pt) for Au-i-Au central collisions at y^S^v = 62.4 GeV 
computed with the global p-p fit [ ^2| ] and the measured reference spectrum [§] are shown, while 
the final Raa{pt) for Au-i-Au central collisions at 62.4 is compared to Raa{pt) at 200 GeV in 
Fig. ^(right) [|23[]. If it weren't already obvious, this should be a lesson to the LHC physicists (and 
management) of the importance of making reference measurements in the same detector for p-p 
collisions at the identical ^/s as the ^/smn of the A-i-A collisions. 

2.1 //^^-suppression — still golden? 

The dramatic difference in suppression of hard-scattering at RHIC compared to SPS fixed 
target cm. energy (^/snn = 17 GeV) stands in stark contrast to suppression, originally thought 
to be the gold-plated signature for deconfinement and the Quark Gluon Plasma (QGP) [^]. Raa 
for suppression is the same, if not identical, at SPS and RHIC (see Fig. ^(left)), thus casting 
a serious doubt on the value of J suppression as a probe of deconfinement. The medium at 
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Figure 7: (left) Comparison of Raa{pt) for K in y/sNN = 62.4 GeV central Au+Au collisions using the 
fit to the previous world 62.4 GeV p-p data or the measured PHENIX reference 62.4 GeV p-p data [^]. 
(right) Final Raa{pt) in central Au+Au collisions at y^s^N — 62.4 and 200 GeV [H. 



RHIC makes 7r"'s nearly vanish but leaves the //*P unchanged compared to lower ^/sNN. One 
possible explanation is that c and c quarks in the QGP recombine to regenerate 7/*P (see Fig. ^ 
(right)), miraculously making the observed Raa equal at SpS and RHIC cm. energies. The good 
news is that such models predict J/^ enhancement (Raa > 1) at LHC energies, which would be 
spectacular, if observed. 
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Figure 8: (left) r'J^ vs centrality (A/pan) at RHIC and SpS energies (right) Predictions for r'J^ in a 
model with regeneration [ 26 1 . 

This leaves us with the interesting question: Will Peter Higgs or Helmut Satz have to wait 
longer at LHC to find out whether they are right? 
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3. The baryon anomaly and xt scaling 



Many RHI physicists tend to treat non-identified cliarged hadrons /j* as if tliey were tt^. 
While this may be a reasonable assumption in p-p collisions, it is clear from Fig. ^-(left) that the 
suppression of non-identified charged hadrons and tt" is very different for 1 < /jj- < 6 GeV/c. 

If the production of high-pr particles in Au-i-Au collisions is the result of hard scattering ac- 
cording to pQCD, then xj scaling should work just as well in Au-i-Au colhsions as in p-p collisions 
and should yield the same value of the exponent nQf[{xT , ^/s) . The only assumption required is that 



the structure and fragmentation functions in Au-i-Au collisions should scale, in which case Eq. |1.2 
still applies, albeit with a G{xt) appropriate for Au-i-Au. In Fig. ^ n^nixT , ^/snn) in Au-i-Au is 
shown for and in peripheral and central collisions, derived by taking the ratio of Ed^o jdp^ 
at a given xj for yjsm =130 and 200 GeV, in each case. [ 27] 
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Figure 9: Power-law exponent neff('^7') for n'^ and /z* spectra in central and peripheral Au+Au collisions at 
= 130 and 200 GeV p7|]. 



The 7r"'s exhibit xj scaling, with the same value of «eff = 6.3 as in p-p collisions, for both 
Au-i-Au peripheral and central collisions. The xj scaling establishes that high-/77 tt" production in 
peripheral and central Au-i-Au collisions follows pQCD as in p-p collisions, with parton distribu- 
tions and fragmentation functions that scale with xj, at least within the experimental sensitivity of 
the data. The fact that the fragmentation functions scale for Tl" in AuH-Au central collisions indi- 
cates that the effective energy loss must scale, i.e. AE{pt) / pr is a constant, which is consistent 
with the constant value of Raa{pt) for pr > 4 GeV/c (Fig. ||-(left)), given that the tt" pr spectrum 
is a pure power-law (Fig. ^(right)-a)). 

The deviation of /j* in Fig. ^ from xj scaling in central Au-i-Au collisions is indicative of and 
consistent with the strong non-scaling modification of particle composition of identified hadrons 
observed in Au-i-Au collisions compared to that of p-p collisions in the range 2.0 <Pt< 4.5 GeV/c, 
where particle production is the result of jet-fragmentation. This is called the Baryon Anomaly. As 
shown in Fig. [TO|-(left) the p/TT^ and p/n^ ratios as a function of pj increase dramatically to 
values ~1 as a function of centrality in Au-i-Au collisions at RHIC [28|. This is nearly an order 



of magnitude larger than had ever been seen previously in either fragmentation of jets in e^e 
collisions or in the average particle composition of the bulk matter in Au-i-Au central collisions [E^] 
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Figure 10: (left) p/n^ and p/n^ as a function of pr and centrality from Au+Au collisions at y^s^N = 
200 GeV compared to other data indicated; (right) Conditional yields, per trigger meson (circles), 
baryon (squares) with 2.5 < pr < 4- GeV/c, of associated mesons with \ .1 < pj < 2.5 GeV/c integrated 
within = ±0.94 radian of the trigger (near side-full) or opposite azimuthal angle (open), for Au+Au 
colUsions at V^aw = 200 GeV @. 



This 'baryon anomaly' was beautifully explained as due to the coalescence of an exponential 
(thermal) distribution of constituent quarks (a.k.a. the QGP) Unfortunately, measurements 
of correlations of in the range 1.7 < pT^ < 2.5 GeV/c associated to identified meson or baryon 
triggers with 2.5 < pt, < 4.0 GeV/c showed the same near side and away side peaks and yields 
(Fig. [I^(right)) characteristic of di-jet production from hard-scattering [32, 30 1, rather than from 
soft coalescence, apparently ruUng out this beautiful model. 

There are still plenty of other models of the baryon anomaly, but none of them are clearly 
definitive. For instance, Stan Brodsky presented at this meeting [33] a higher twist model of the 
baryon anomaly as the result of the reaction q + q P + q- This predicts an isolated proton with 
no same-side jet, but with an opposite jet, a clear and crucial test. Another test (from the CIM [^) 
is that Weff 8 for these protons. This effect will be emphasized in central collisions because the 
higher twist subprocesses have 'small size' and are 'color transparent' so they propagate through 
the nuclear medium without absorption. This is consistent with the reduced near-side correlation 
to baryon triggers compared to meson triggers shown for the most central collisions in Fig. IC- 
(right) [ 30 1 ; but definitive detection of isolated p or pin central A+A collisions or precision mea- 
surements of XT scaling for p and as a function of centrality remain very interesting projects for 
the future. 



Guy Paic recently [ ]34| ] claimed to explain the baryon anomaly by simple radial flow, which 
occurs late in the expansion, even in the hadronic phase. The radial flow velocity Lorentz-boosts 
the heaver protons to larger pT than the lighter pions. Also, protons have a shorter formation time 
than pions so they may participate in the radial flow even if they result from parton fragmentation. 
If Guy is correct, then it is time to re-examine the subject of the elliptic flow (V2) and pj spectra 
of identified hadrons which was very popular several years ago. The last time I looked, in the 
PHENIX White Paper [35], the hydro models with radial flow could either explain the pj spectra 
of p and n or the V2 but not both. If the steadily improving models can now explain both V2 and the 
Pt spectra, this would spell the end of the 'baryon anomaly'. However, the anomalously large p/n 
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ratio would remain a signature property of the inclusive identified particle pr spectra in central 
Au+Au collisions. 

In ttiis vein, I was temporarily thrown for a loop by a result from a STAR presentation at 
Quark Matter 2008, as reported by Marco van Leeuwen at Hard Probes 2008 [ ^ ] (and Christine 
Nattrass [37] at this meeting), which measured the particle composition in the near-side jet and 
'ridge' and seemed to indicate that the p/n (or baryon/meson) ratio in the near-side jet was not 
anomalous. I first thought that this disagreed with everything that I said previously in this section 
about the 'baryon anomaly'. To quote Marco [36], "the p/n ratio in the ridge is similar to the 
inclusive p/n ratio in Au-i- Au events, which is much larger than in p-i-p events. The p/n ratio in 
the jet-like peak is similar to the inclusive ratio in p -i- p events." However, at a Ridge workshop 
at BNL [38], I found out that what STAR really meant to say was that "The p/n ratio of the 
conditional yield for near side-correlations associated to an h"^ trigger with pr, > 4 GeV/c is similar 
in the jet-like peak to the inclusive ratio in p -i- p events; while the p/n ratio in the ridge is similar 
to the inclusive p/n ratio in Au-i- Au events, which is much larger than in p-np events." — i.e. STAR 
was talking about associated yields to the trigger /z^ and did not include the triggering particle in 
the yield. Hence there is no disagreement. In fact, the STAR result is actually in agreement with a 
recent PHENIX measurement [39] of the ratio of the associated baryon and meson conditional near- 
side yields from an /j^ trigger with 2.5 < pr, < 4 GeV/c in Au-i-Au collisions at ^/sf^ = 200 GeV. 

To reiterate, the anomalously large p/n ratio remains a signature property of the inclusive 
identified particle pr spectra in central Au-i-Au collisions. In fact, Christine Nattrass' [|^] thorough 
demonstration at this meeting of the properties and particle composition of 'the ridge' convinced 
me that that the ridge is nothing other than a random coincidence of any trigger particle and the 
background or bulk of inclusive particles, which, as Stan Brodsky commented, is generally biased 
towards the trigger direction due to the effect. Naturally, several details such as the azimuthal 
width of the ridge still remain to be explained. In my opinion, a key test of this idea is that a 
ridge of same side correlations with large Arj should exist for direct-y triggers, or, in the PHENIX 
acceptance, the same-side correlation to a direct-/ should exist at the same rate and azimuthal width 
as the ridge we observed in tt" or inclusive 7 same-side correlations. [ pO| ] 



4. Direct photons at RHIC — Thermal photons? 

4.1 Internal Conversions — the first measurement anywhere of direct photons at low pr 

Internal conversion of a photon from tt" and rj decay is well-known and is called Dalitz de- 
cay [pr]]. Perhaps less well known in the RHI community is the fact that for any reaction (e.g. 
q + g ^ Y + q) in which a real photon can be emitted, a virtual photon (e.g. e^e^ pair of mass 
rriee > 2me) can also be emitted. This is called internal-conversion and is generally given by the 
Kroll-Wada formula [||, ||]: 



1 dNee 
Nj dniee 




'""'1 + ^) , (4.1) 
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where M is the mass of the decaying meson or the effective mass of the emitting system. The domi- 
nant terms are on the first line of Eq. |4.l[ the characteristic 1 /m^g dependence; and the cutoff of the 
spectrum for > M (Fig. 0-(Ieft)) [Q. Since the main background for direct-single- 7 produc- 
tion is a photon from tt" ^ 7 + 7, selecting mgg> 100 MeV effectively reduces the background by 
an order of magnitude by eliminating the background from Dalitz decay, ^ j + e^ + e^, at 
the expense of a factor ~ 1000 in rate. This allows the direct photon measurements to be extended 
(for the first time in both p-p and Au-i-Au collisions) below the value of ~ 4 GeV/c, possible with 
real photons, down to pj = I GeV/c (Fig. [ri]-(right)) [^], which is a real achievement. The solid 
lines on the p-p data are QCD calculations which work down to pr =2 GeV/c. The dashed line is 
a fit of the p-p data to the modified power law B{1 + pj/b)~", used in the related Drell-Yan [44] 
reaction, which flattens as pj — > 0. 

The relatively flat, non-exponential, spectra for the direct-7 and Drell-Yan reactions as pj —>^0 
is due to the fact that there is no soft-physics production process for them, only production via 
the partonic subprocesses, g + q —>^ 7+1 and q + q + e^, respectively. This is quite distinct 
from the case for hadron production, e.g. tt", where the spectra aie exponential as — > in 
p-p collisions (Fig. [l]) due to soft-production processes, as well as in Au-i-Au collisions. Thus, 
for direct-7 iri Au-i-Au collisions, the exponential spectrum of excess photons above the {Taa) 
extrapolated p-p fit is unique and therefore suggestive of a thermal source. 
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Figure 11: (left) Invariant mass (m^^^-) distribution of e+e pairs from Au+Au minimum bias events for 



1.0 < pt < 1.5 GeV/c [[4-3p. Dashed lines are Eq. 4-.1 for the mesons indicated. Blue solid line is fc{m), the 
total di-electron yield from the sum of contributions or 'cocktail' of meson Dalitz decays; Red solid line is 
fdir{fn) the internal conversion m^+g- spectrum from a direct-photon (M >> mg+g ). Black solid line is a fit 
of the data to the sum of cocktail plus direct contributions in the range 80 < m,+,- < 300 MeV/c2. (right) 
Invariant cross section (p-p) or invariant yield (Au+Au) of direct photons as a function of pj Filled 
points are from virtual photons, open points from real photons. 
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4.2 Low pt vs high pr direct- 7 — Learn a lot from a busy plot 



The unique behavior of direct-7 at low pr in Au+Au relative to p+p compared to any other par- 
ticle is more dramatically illustrated by examining the Raa of all particles measured by PHENIX in 
central Au+Au collisions at 



= 200 GeV (Fig. |12|) For the entire region pr < 20 GeV/c 
so far measured at RHIC, apart from the p + p which are enhanced in the region 2 < pj <4 GeV/c 
('the baryon anomaly'), the production of no other particle is enhanced over point-like scaling. 
The behavior of Raa of the low pj < 2 GeV/c direct-/ is totally and dramatically different from 
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Figure 12: Nuclear Modification Factor, Raa{pt) for all identified particles so far measured by PHENIX in 
central Au+Au collisions at ^/snn — 200GeV. [|7[ 

all the other particles, exhibiting an order of magnitude exponential enhancement as pr — > 0. This 
exponential enhancement is certainly suggestive of a new production mechanism in central Au+Au 
collisions different from the conventional soft and hard particle production processes in p-p colli- 



sions and its unique behavior is attributed to thermal photon production by many authors. [48] 



4.2.1 Direct photons and mesons up to pj = 20 GeV/c 

Other instructive observations can be gleaned from Fig. The tt" and 77 continue to track 
each other to the highest pj. At lower px, the <p meson tracks the A"^ very well, but with a different 
value of Raa{pt) than the n^, while at higher pT,the (p and (0 vector mesons appear to track each 
other. Interestingly, the 7/*P seems to track the tt" for < /jj- < 4 GeV/c; and it will be interesting 
to see whether this trend continues at higher pj. 
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The direct-7's also show something interesting at high pr which might possibly indicate trou- 
ble ahead at the LHC. With admittedly large systematic errors, which should not be ignored, the 
direct-7 appear to become suppressed for pj > 14 GeV/c with a trend towards equality with 

y 

for pj ~ 20 GeV. Should become equal to 7?^, it would imply that the energy loss in the final 
state is no longer a significant effect for pj>20 GeV/c and that the equal suppression of direct-/ 
and is due to the initial state structure functions. If this were true, it could mean that going 
to much higher pr would not be useful for measurements of parton suppression. In this vein, 
the new EPS09 structure functions for quarks and gluons in nuclei were presented at this meet- 
ing [|49[], which represented the best estimate of shadowing derived by fitting all the DIS data in 
IJ.{e) — A scattering as well as including in the fit, notably, the PHENIX data in d-i-Au and p-p 
as a function of centrality. Clearly, improved measurements of both direct-7 and tt" in the range 
10 < Pr < 20 GeV/c are of the utmost importance for both the RHIC and LHC programs. 



5. Precision measurements, key to the next step in understanding 

There are many different models of parton suppression with totally different assumptions 
which all give results in agreement with the PHENIX measurement R^ 0.20 for 4 < pj < 
20 GeV/c in Au-i-Au central collisions. In PHENIX, Jamie Nagle got all theorists to send us pre- 
dictions as a function of their main single parameter that characterizes the medium in order to do 
precision fits to the latest PHENIX data including the correct treatment of correlated experi- 
mental systematic errors (Fig. 13 ) [|50|]. Systematic uncertainties of the theory predictions were 




Figure 13: a) (left) PHENIX k° Raa{pt) for Au+Au central (0-5%) collisions at ./^n = 200 [50| compared 



to PQM model predictions [51 1 as a function of {q). The thick red line is the best fit. b) (center) Values of 



Raa at pt = 20 GeV/c as a function of {q) in the PQM model [|T| corresponding to the lines the left panel, 
c) (right) same as b) but on a log-log scale, with fit. 



not considered. 

The large value of the transport coefficient (q 



13.2+^;^ GeV^/fm from the best fit 



to the PQM model [51] (where /X is the average 4-momentum transfer to the medium per mean 
free path A) is a subject of some debate in both the more fundamental QCD community [|5^] and 



the more phenomenological community [53|. For instance it was stated in Ref. [53] that "the 



dependence of Raa on q becomes weaker as q increases" as is clear from Fig. 13 b. It was also 
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asserted that "when the values of the time-averaged transport coefficient q exceeds 5 GeV^/fm, 
Raa gradually loses its sensitivity." That statement also appeared reasonable. However, given the 
opportunity of looking at a whole range of theoretical predictions (kindly provided by the PQM 
authors [|5T|]) rather than just the one that happens to fit the data, we experimentalists learned 
something about the theory that was different from what the theorists emphasized. By simply 
looking at the PQM predictions on a log-log plot (Fig. |l3b), it became evident that the PQM 



prediction could be parameterized as Raa [pt = 20 GeV/c] = 0.75 / \/q{ GeV^/fm) over the range 
5 < ^ < 100 GeV^/fm. This means that in this range, the fractional sensitivity to q is simply 
proportional to the fractional uncertainty in Raa, i-S- ^q/q = 2.0 x ARaa/Raa, so that improving 
the precision of Raa e.g. in the range 10 < < 20 GeV/c will lead to improved precision on (17). 
This is a strong incentive for experimentalists. Similarly, this should give the theorists incentive to 
improve their (generally unstated) systematic uncertainties. 

5.1 Raa vs. the reaction plane 

Another good synergy between experimentalists and theorists is the study of Raa as a function 
of angle to the reaction plane and centrality in order to understand the effect of varying the initial 
conditions (centrality) and the path length through the medium (angle). When PHENIX first pre- 



sented results on Raa{pt) vs. the angle A0 to the reaction plane [54] there was a reaction from 
the flow community that this is nothing other than a different way to present the anisotropic flow, 
V2. This is strictly not true for two reasons: 1) V2 measurements are relative while Raa{^^,Pt) 
is an absolute measurement including efficiency, acceptance and all other such corrections; 2) if 
and only if the angular distribution of high pr suppression around the reaction plane were sim- 
ply a second harmonic so that all the harmonics other than V2 vanish (and why should that be?) 
then Raa{^<^ -iPt) / Raa{pt) = 1 +2v2Cos2A0. Nevertheless, whatever the actual form of the an- 
gular distribution, it is true that Raa{^^ , Pt) / Raa{pt) = dN{/S.^,pT)/d/S.^ / {dN (Acj) , pr) / dA(j))) 
but without the absolute values it is impossible to tell whether Raa{A<P,pt) approaches or exceeds 
1 (or any other value) at some value of A(j). 

For instance in a new result this year, PHENIX has observed a striking difference in the be- 
havior of the dependence of the in-plane Raa{A<P ^ Q,pt) for tt" as a function of centrality N^an 
compared to the dependence of the Raa{A<^ ~ 7i/2,pj) in the direction perpendicular to the reac- 



tion plane | |55p (Fig. |14|). The Raa perpendicular to the reaction plane is relatively constant with 
centrality, while the predominant variation of Raa with centrality comes in the direction parallel to 
the reaction plane which shows a strong centrality dependence. This is a clear demonstration of the 
sensitivity of Raa to the length traversed in the medium which is relatively constant as a function of 
centrality perpendicular to the reaction plane but depends strongly on the centrality parallel to the 
reaction plane. This is a fantastic but reasonable result and suggests that tests of the many models 
of energy loss should concentrate on comparing the centrality dependence for directions parallel to 
the reaction plane, where the length traversed depends strongly on centrality, compared to perpen- 
dicular to the reaction plane, where the length doesn't change much with centrality, before tackling 
the entire angular distribution. 

The theorists have not been idle on this issue and are making great strides by attempting to put 
all the theoretical models of jet quenching into a common nuclear geometrical and medium evo- 
lution formalism so as to get an idea of the fundamental differences in the models [|5^ "evaluated 
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Figure 14: Nuclear Modification Factor, in reaction-plane bins as a function of pr and centrality (A^part) 
in Au+Au collisions at ^/snn = 200 GeV |^]. Filled circles represent ^^4(0 < A(p < 15°) (in-plane), open 
squares RaaOS < A(j) < 90°) (out-of -plane) and filled triangles /?,u(30 < A<j) < 45°). 



on identical media, initial state and final fragmentation. Tlie only difference in models will be in 
the Eloss kernel.". The different models all agreed with the measured Raa{pt) (Fig ]_5i); but 
the agreement with the measured Raa{^^,Pt) as shown by the 7?^ (out )//?aa (in) ratio is not very 
good (Fig 15b). Hopefully the latest PHENIX results [55] (Fig. [T^ will suggest the way for further 
improvement. 
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Figure 15: a) (left) Three theoretical models [|6|] compared to PHENIX R'^{pt) in AuH-Au collisions at 
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6. Do direct-e^ from Heavy Flavors indicate one or two theoretical crises? 

PHENIX was specifically designed to be able to detect charm particles via direct single-e* 
from their semileptonic decay. Fig. 16 a shows our direct single-e* measurement in p-p collisions 

at ^/s-■ 



200 GeV pTp in agreement with a QCD calculation [ |58| ] of c and b quarks as the source 
of the direct single-e^ (heavy-flavor e^). In Au+Au collisions, a totally unexpected result was ob- 




Figure 16: a) (left) Invariant cross section of direct single-e* in p-p collisions ^% compared to theoretical 
predictions from c and b quark semileptonic decay. [|8]| b) (right) Raa as a function of pr for direct single- 
e"*^ nP and r\ in Auh-Au central (0-10%) collisions at = 200 GeV. 



served. The direct single-e^ from heavy quarks are suppressed the same as the and x] from light 
quarks (and gluons) in the range 4 < pT <9 GeV/c where b and c contributions are roughly equal 
(Fig. |I^) ||59|]. This strongly disfavors the QCD energy-loss explanation of jet-quenching because, 
naively, heavy quarks should radiate much less than light quarks and gluons in the medium; but 
opens up a whole range of new possibilities including string theory []60|]. 

6.1 Zichichi to the rescue? 

In September 2007, I read an article by Nino Zichichi, "Yukawa's gold mine", in the CERN 
Courier, taken from his talk at the 2007 International Nuclear Physics meeting in Tokyo, Japan, 
in which he proposed: "the reason why the top quark appears to be so heavy (around 200 GeV) 
could be the result of some, so far unknown, condition related to the fact that the final state must be 
QCD-colourless. We know that confinement produces masses of the order of a giga-electron-volt. 
Therefore, according to our present understanding, the QCD colourless condition cannot explain 
the heavy quark mass. However, since the origin of the quark masses is still not known, it cannot 
be excluded that in a QCD coloured world, the six quarks are all nearly massless and that the 
colourless condition is 'flavour' dependent." 

Nino's idea really excited me, even though, or perhaps, because, it appeared to overturn two 
of the major tenets of the Standard Model since it seemed to imply that: QCD isn't flavor blind; the 
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masses of quarks aren't given by the Higgs mechanism. Massless b and c quarks in a color-charged 
medium would be the simplest way to explain the apparent equality of gluon, light quark and heavy 
quark suppression indicated by the equality of Raa for and Raa of direct single-e^ in regions 
where both c and b quarks dominate. Furthermore RHIC and LHC-Ions are the only place in the 
Universe to test this idea. Nino's idea seems much more reasonable to me than the string theory 
explanations of heavy-quark suppression (especially since they can't explain light-quark suppres- 
sion). Nevertheless, just to be safe, I asked some distinguished theorists what they thought, with 
these results: "Oh, you mean the Higgs field can't penetrate the QGP" (Stan Brodsky); "You mean 
that the propagation of heavy and light quarks through the medium is the same" (Rob Pisarski); 
"The Higgs coupling to vector bosons 7, W, Z is specified in the standard model and is a funda- 
mental issue. One big question to be answered by the LHC is whether the Higgs gives mass to 
fermions or only to gauge bosons. The Yukawa couplings to fermions are put in by hand and are 
not required" "What sets fermion masses, mixings?" (Chris Quigg-Moriond2008); "No change in 
the f-quark, W, Higgs mass relationship if there is no Yukawa coupling: but there could be other 
changes" (Bill Marciano). 

Nino proposed to test his idea by shooting a proton beam through a QGP formed in a Pb-i-Pb 
collision at the LHC and seeing the proton 'dissolved' by the QGP. My idea is to use the new 
PHENIX vertex detector, to be installed in 2010, to map out, on an event-by-event basis, the di- 
hadron correlations from identified b,b di-jets, identified c,c di-jets, which do not originate from 
the vertex, and light quark and gluon di-jets, which originate from the vertex and can be measured 
with TT^-hadron correlations. These measurements will confirm in detail (or falsify) whether the 
different flavors of quarks behave as if they have the same mass in a color-charged medium. De- 
pending when the LHC-Ions starts, it is conceivable that ALICE or another LHC experiment with 
a good vertex detector could beat RHIC to the punch, since this measurement compares the energy 
loss of light and heavy quarks and may not need p-p comparison data. 

If Nino's proposed effect is true, that the masses of fermions are not given by the Higgs, and 
we can confirm the effect at RHIC or LHC-Ions, this would be a case where Relativistic Heavy Ion 
Physics may have something unique to contribute at the most fundamental level to the Standard 
Model — a "transformational discovery." Of course the LHC or Tevatron could falsify this idea by 
finding the Higgs decay to b,b at the expected rate in p-p collisions. 



7. Soft physics projections for LHC 



Some soft physics issues at LHC are also very interesting to me. Marek Gazdzicki has popular- 



ized 3 features from the NA49 results [ pip at the CERN SpS fixed target heavy ion program, which 
he calls 'the kink', 'the horn' and 'the step'. I believe that 'the kink' is certainly correct (Fig. |l7|a) 
and has relevance to the LHC program. The 'kink' reflects the fact that the wounded nucleon model 

TV ~ 20 GeV where it was discovered [|64l, |65ll and fails above and 



(WNM) [63] works only at 



below y/sNN ~ 20 GeV: wounded projectile nucleons below 20 GeV at mid-rapidity 
projectile quarks (AQM) 



Tp, 31 GeV and above |68, 69] 



6p; wounded 

This led me to speculate that maybe the 
charged-particle multiplicity or sum-transverse energy might be the only quantity to exhibit point- 
like Ncow scaling at LHC energies. However, Wit Busza's prediction for the charged multiplicity 
per participant pair increased by the same ratio from RHIC to LHC in both A-i-A and p-p collisions. 
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Figure 17: a) (left) Gazdzciki's plot [ ]6l| ] of pions/participant vs. F = {^/snn — Int^)^^^ / ^/snn^^"^ 



/Snn 



1/2 



in A+A and p+p collisions, b) (right) Wit Busza's prediction for the number of charged parti- 



cles per participant pair vs. In y/s(GeW) [p2t. 



which implies that the AQM will still work at the LHC. This makes me think that A/^^oii scaling for 
soft-processes at LHC is unlikely. 

A more interesting soft physics issue for the LHC concerns the possible increase of the an- 
isotropic flow V2 beyond the 'hydrodynamic limit'. Wit Busza's extrapolation [ |62| ] of V2 to the LHC 
energy is shown in Fig. |l^a, a factor of L6 increase from RHIC. A previous paper by NA49 [7C] 
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Figure 18: a) (left) Busza's extrapolation of V2 to LHC |^]. b) (center) V2/E vs 'Bjorken multiplicity 
density', {I /S)dNch/dy [f70|]. c) (right) 'Hydro Limit' calculated in viscous Hydrodynamics for several 
values of the initial energy density eo 0] ■ 



which compared V2 measurements from AGS and CERN fixed target experiments to RHIC as a 
function of the 'Bjorken multiplicity density', dncb/drj /S, where 5 = is the overlap area of the 
collision zone, showed an increase in V2/£ from fixed target energies to RHIC leading to a "hydro 
limit", where e is the eccentricity of the collision zone (Fig. [T8|b). This Umit was confirmed in a 
recent calculation using viscous relativistic hydrodynamics [ f7T| ] which showed a clear hydro-limit 
of V2/£ = 0.20 (Fig. 18 c). This limit is sensitive to the ratio of the viscosity /entropy density, the 
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now famous rj/s, but negligibly sensitive to the maximum energy density of the collision, so I 
assume that this calculation would give a hydro-limit at the LHC not too different from RHIC, 
V2/£ ~ 0.20. Busza's extrapolation of a factor of 1.6 increase in V2 from RHIC to LHC combined 
with V2/£ from Fig. |^ gives V2/e = 0.32 at LHC. In my opinion this is a measurement which can 
be done to high precision on the first day of Pb+Pb collisions at the LHC, since it is high rate and 
needs no p-p comparison data. Personally, I wonder what the hydro aficionados would say if both 
Heinz and Busza's predictions were correct? 
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